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Abstract 

I review the current state of determinations of the Hubble constant, 
which gives the length scale of the Universe by relating the expansion ve- 
locity of objects to their distance. In the last 20 years, much progress has 
been made and estimates now range between 60 and 75 km s _1 Mpc -1 , 
with most now between 70 and 75 km s _1 Mpc -1 , a huge improvement 
over the factor-of-2 uncertainty which used to prevail. Further improve- 
ments which gave a generally agreed margin of error of a few percent 
rather than the current 10% would be vital input to much other interest- 
ing cosmology. There are several programmes which are likely to lead us 
to this point in the next 10 years. 
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1 Introduction 



1.1 A brief history 

The last century saw an expansion in our view of the world from a static, 
Galaxy-sized Universe, whose constituents were stars and "nebulae" of unknown 
but possibly stellar origin, to the view that the observable Universe is in a state 
of expansion from an initial singularity over ten billion years ago, and contains 
approximately 100 billion galaxies. This paradigm shift was summarised in a 
famous debate between Shapley and Curtis in 1920; summaries of the views of 
each protagonist can be found in [57] and [140] . 

The historical background to this change in world view has been extensively 
discussed and whole books have been devoted to the subject of distance mea- 
surement in astronomy [125] . At the heart of the change was the conclusive 
proof that what we now know as external galaxies lay at huge distances, much 
greater than those between objects in our own Galaxy. The earliest such dis- 
tance determinations included those of the galaxies NGC 6822 |61| . M33 [B5] 
and M31 E>3|. 

As well as determining distances, Hubble also considered redshifts of spectral 
lines in galaxy spectra which had previously been measured by Slipher in a series 
of papers 142, 143 . If a spectral line of emitted wavelength Ao is observed at 
a wavelength A, the redshift z is defined as 

z = A/A - 1. (1) 

For nearby objects, the redshift corresponds to a recession velocity v which for 
nearby objects is given by a simple Doppler formula, v = cz. Hubble showed 
that a relation existed between distance and redshift (see FigureQ]); more distant 
galaxies recede faster, an observation which can naturally be explained if the 
Universe as a whole is expanding. The relation between the recession velocity 
and distance is linear, as it must be if the same dependence is to be observed 
from any other galaxy as it is from our own Galaxy (see Figure [2]) . The pro- 
portionality constant is the Hubble constant Hq, where the subscript indicates 
a value as measured now. Unless the Universe's expansion does not accelerate 
or decelerate, the slope of the velocity-distance relation is different for observers 
at different epochs of the Universe. 

Recession velocities are very easy to measure; all we need is an object with an 
emission line and a spectrograph. Distances are very difficult. This is because 
in order to measure a distance, we need a standard candle (an object whose 
luminosity is known) or a standard ruler (an object whose length is known), 
and we then use apparent brightness or angular size to work out the distance. 
Good standard candles and standard rulers are in short supply because most 
such objects require that we understand their astrophysics well enough to work 
out what their luminosity or size actually is. Neither stars nor galaxies by 
themselves remotely approach the uniformity needed; even when selected by 
other, easily measurable properties such as colour, they range over orders of 
magnitude in luminosity and size for reasons that are astrophysically interesting 
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Figure 1: Hubble's original diagram of distance to nearby galaxies, derived 
from measurements using Cepheid variables, against velocity, derived from red- 
shift \b °J§ . The Hubble constant is the slope of this relation, and in this diagram 
is a factor of nearly 10 steeper than currently accepted values. 




Figure 2: Illustration of the Hubble law. Galaxies at all points of the square grid 
are receding from the black galaxy at the centre, with velocities proportional to 
their distance away from it. From the point of view of the second, green, galaxy 
two grid points to the left, all velocities are modified by vector addition of its 
velocity relative to the black galaxy (red arrows). When this is done, velocities 
of galaxies as seen by the second galaxy are indicated by green arrows; they all 
appear to recede from this galaxy, again with a Hubble-law linear dependence of 
velocity on distance. 
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but frustrating for distance measurement. The ideal Hq object, in fact, is one 
which involves as little astrophysics as possible. 

Hubble originally used a class of stars known as Cepheid variables for his 
distance determinations. These are giant blue stars, the best known of which is 
aUMa, or Polaris. In most normal stars, a self-regulating mechanism exists in 
which any tendency for the star to expand or contract is quickly damped out. 
In a small range of temperature on the Hertzsprung-Russell (H-R) diagram, 
around 7000-8000 K, particularly at high luminosity, this does not happen 
and pulsations occur. These pulsations, the defining property of Cepheids, have 
a characteristic form, a steep rise followed by a gradual fall, and a period which 
is directly proportional to luminosity. The period-luminosity relationship was 
discovered by Leavitt [86] by studying a sample of Cepheid variables in the Large 
Magellanic Cloud (LMC). Because these stars were known to be all at the same 
distance, their correlation of apparent magnitude with period therefore implied 
the P-L relationship. 

The Hubble constant was originally measured as 500 km s _1 Mpc -1 [53] 
and its subsequent history was a more-or-less uniform revision downwards. In 
the early days this was caused by biat@ in the original samples [8j, confusion 
between bright stars and Hn regions in the original samples [55] 1131] and dif- 
ferences between type I and II Cepheidfd 4J. In the second half of the last 
century, the subject was dominated by a lengthy dispute between investiga- 
tors favouring values around 50 km s _1 Mpc -1 and those preferring higher val- 
ues of 100 km s -1 Mpc -1 . Most astronomers would now bet large amounts 
of money on the true value lying between these extremes, and this review is 
an attempt to explain why and also to try and evaluate the evidence for the 
best-guess (2007) current value. It is not an attempt to review the global his- 
tory of Hq determinations, as this has been done many times, often by the 
original protagonists or their close collaborators. For an overall review of this 

1 This is known in the literature as the "instability strip" and is almost, but not quite, 
parallel to the luminosity axis on the H-R diagram; brighter Cepheids have slightly lower 
temperatures. The instability strip has a finite width, which causes a small degree of dispersion 
in period— luminosity correlations among Cepheids. 

2 There are numerous subtle and less-subtle biases in distance measurement; see 11511 for 
a blow-by-blow account. The simplest bias, the "classical" Malmquist bias, arises because, 
in any population of objects with a distribution in intrinsic luminosity, only the brighter 
members of the population will be seen at large distances. The result is that the inferred 
average luminosity is greater than the true luminosity, biasing distance measurements towards 
the systematically short. The Behr bias [8] from 1951 is a distance-dependent version of the 
Malmquist bias, namely that at higher distances, increasingly bright galaxies will be missing 
from samples. This leads to an overestimate of the average brightness of the standard candle 
which becomes worse at higher distance. 

3 Cepheids come in two flavours: type I and type II, corresponding to population I and II 
stars. Population II stars are the first generation of stars, which formed before the enrichment 
of the ISM by detritus from earlier stars, and Population I stars like the Sun are the later 
generation which contain significant amounts of elements other than hydrogen and helium. 
The name "Cepheid" derives from the fact that the star <5 Cephei was the first to be identified 
(by Goodricke in 1784). Population II Cepheids are sometimes known as W Virginis stars, 
after their prototype, W Vir, and a W Vir star is typically a factor of 3 fainter than a classical 
Cepheid of the same period. 



4 



process see, for example, |161] and |149| ; see also data compilations and re- 
views by Huchra (http ://cfa-www. harvard. edu^huchra/hubble) and Allen 
( jhttp : //www, i nstitute-of -brilliant-failure s . com/p . 

1.2 A little cosmology 

The expanding Universe is a consequence, although not the only possible conse- 
quence, of general relativity coupled with the assumption that space is homoge- 
neous (that is, it has the same average density of matter at all points at a given 
time) and isotropic (the same in all directions). In 1922 Friedman [47] showed 
that given that assumption, we can use the Einstein field equations of general 
relativity to write down the dynamics of the Universe using the following two 
equations, now known as the Friedman equations: 

kc 2 , (2) 
inG(p + 3p/c 2 ) + ±A. (3) 

Here a = a(t) is the scale factor of the Universe. It is fundamentally related to 
rcdshift, because the quantity (1 + z) is the ratio of the scale of the Universe 
now to the scale of the Universe at the time of emission of the light (ao/a). 
A is the cosmological constant, which appears in the field equation of general 
relativity as an extra term. It corresponds to a universal repulsion and was 
originally introduced by Einstein to coerce the Universe into being static. On 
Hubble's discovery of the expansion of the Universe, he removed it, only for it to 
reappear seventy years later as a result of new data [1161 1123] (see also [21] for 
a review). A: is a curvature term, and is —1, 0, or +1, according to whether the 
global geometry of the Universe is negatively curved, spatially flat, or positively 
curved, p is the density of the contents of the Universe, p is the pressure and 
dots represent time derivatives. For any particular component of the Universe, 
we need to specify an equation for the relation of pressure to density to solve 
these equations; for most components of interest such an equation is of the form 
p = wp. Component densities vary with scale factor a as the Universe expands, 
and hence vary with time. 

At any given time, we can define a Hubble parameter 

H[t) = a/a, (4) 

which is obviously related to the Hubble constant, because it is the ratio of an 
increase in scale factor to the scale factor itself. In fact, the Hubble constant 
Hq is just the value of H at the current time. 

If A = 0, we can derive the kinematics of the Universe quite simply from the 
first Friedman equation. For a spatially flat Universe k = 0, and we therefore 
have 

P = Pc -^ (5) 



a 2 - -{^Gp + K)a 2 = - 

a 
a 
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where p c is known as the critical density. For Universes whose densities are 
less than this critical density, k < and space is negatively curved. For such 
Universes it is easy to see from the first Friedman equation that we require a > 0, 
and therefore the Universe must carry on expanding for ever. For positively 
curved Universes (k > 0), the right hand side is negative, and we reach a 
point at which a = 0. At this point the expansion will stop and thereafter go 
into reverse, leading eventually to a Big Crunch as a becomes larger and more 
negative. 

For the global history of the Universe in models with a cosmological constant, 
however, we need to consider the A term as providing an effective acceleration. 
If the cosmological constant is positive, the Universe is almost bound to expand 
forever, unless the matter density is very much greater than the energy density 
in cosmological constant and can collapse the Universe before the acceleration 
takes over. (A negative cosmological constant will always cause recollapse, but 
is not part of any currently likely world model) . [21] provides further discussion 
of this point. 

We can also introduce some dimensionless symbols for energy densities in the 
cosmological constant at the current time, J7a = A/(3i?g), and in "curvature 
energy", 0,% = kc 2 /Hq. By rearranging the first Friedman equation we obtain 

^! = £_ _ n k a- 2 + fi A . (6) 

H Pc 

The density in a particular component of the Universe A, as a fraction of 
critical density, can be written as 

Px/Pc = V x a a , (7) 

where the exponent a represents the dilution of the component as the Uni- 
verse expands. It is related to the w parameter defined earlier by the equation 
a = —3(1 + w); Equation Q holds provided that w is constant. For ordinary 
matter a = —3, and for radiation a — —4, because in addition to geometrical 
dilution the energy of radiation decreases as the wavelength increases, in ad- 
dition to dilution due to the universal expansion. The cosmological constant 
energy density remains the same no matter how the size of the Universe in- 
creases, hence for a cosmological constant we have a = and w = — 1. w = — 1 
is not the only possibility for producing acceleration, however; any general class 
of "quintessence" models for which w < — | will do. Moreover, there is no 
reason why w has to be constant with redshift, and future observations may be 
able to constrain models of the form w — Wo + w\z. The term "dark energy" is 
usually used as a general description of all such models, including the cosmolog- 
ical constant; in most current models, the dark energy will become increasingly 
dominant in the dynamics of the Universe as it expands. 
In the simple case, 

Y,n x + n A + n k = i (8) 

x 
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by definition, because Slfe = implies a flat Universe in which the total energy 
density in matter together with the cosmological constant is equal to the critical 
density. Universes for which Qk is almost zero tend to evolve away from this 
point, so the observed near-flatness is a puzzle known as the "flatness problem" ; 
the hypothesis of a period of rapid expansion known as inflation in the early 
history of the Universe predicts this near-flatness naturally. 

We finally obtain an equation for the variation of the Hubble parameter with 
time in terms of the Hubble constant (see e.g. [HI]), 

H 2 = H 2 (Q A + Q m a~ 3 + fl r a- 4 - f^a" 2 ), (9) 

where Q r represents the energy density in radiation and Q m the energy density 
in matter. 

We can define a number of distances in cosmology. The most important 
for present purposes are the angular diameter distance Z?a, which relates the 
apparent angular size of an object to its proper size, and the luminosity dis- 
tance -Dl = (1 + z) 2 Da, which relates the observed flux of an object to its 
intrinsic luminosity. For currently popular models, the angular diameter dis- 
tance increases to a maximum as z increases to a value of order 1, and decreases 
thereafter. Formulae for, and fuller explanations of, both distances are given 
by [56]. 
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2 Local Methods and Cepheid Variables 



2.1 Preliminary remarks 

As we have seen, in principle a single object whose spectrum reveals its re- 
cession velocity, and whose distance or luminosity is accurately known, gives 
a measurement of the Hubble constant. In practice, the object must be far 
enough away for the dominant contribution to the motion to be the velocity 
associated with the general expansion of the Universe (the "Hubble flow"), as 
this expansion velocity increases linearly with distance whereas other nuisance 
velocities, arising from gravitational interaction with nearby matter, do not. For 
nearby galaxies, motions associated with the potential of the local environment 
are about 200-300 km s _1 , requiring us to measure distances corresponding 
to recession velocities of a few thousand km s" 1 or greater. These recession 
velocities correspond to distances of at least a few tens of Mpc. 

For large distances, corresponding to redshifts approaching 1, the relation 
between spectroscopically measured redshift and luminosity (or angular diame- 
ter) distance is no longer linear and depends on the matter density f2 m and dark 
energy density as well as the Hubble constant. This is less of a problem, 
because as we shall see in Section [31 these parameters are probably at least as 
well determined as the Hubble constant itself. 

Unfortunately, there is no object, or class of object, whose luminosity can 
be determined unambiguously in a single step and which can also be observed 
at distances of tens of Mpc. The approach, used since the original papers by 
Hubble, has therefore been to measure distances of nearby objects and use 
this knowledge to calibrate the brightness of more distant objects compared 
to the nearby ones. This process must be repeated several times in order to 
bootstrap one's way out to tens of Mpc, and has been the subject of many 
reviews and books (see e.g. |125j ). The process has a long and tortuous history, 
with many controversies and false turnings, and which as a by-product included 
the discovery of a large amount of stellar astrophysics. The astrophysical content 
of the method is a disadvantage, because errors in our understanding propagate 
directly into errors in the distance scale and consequently the Hubble constant. 
The number of steps involved is also a disadvantage, as it allows opportunities 
for both random and systematic errors to creep into the measurement. It is 
probably fair to say that some of these errors are still not universally agreed 
on. The range of recent estimates is from 60 to 75 km s _1 Mpc -1 , and the 
reasons for the disagreements (in many cases by different analysis of essentially 
the same data) are often quite complex. 

2.2 Basic principle 

We first outline the method briefly, before discussing each stage in more detail. 
Nearby stars have a reliable distance measurement in the form of the parallax 
effect. This effect arises because the earth's motion around the sun produces 
an apparent shift in the position of nearby stars compared to background stars 



8 



at much greater distances. The shift has a period of a year, and an angular 
amplitude on the sky of the Earth-Sun distance divided by the distance to the 
star. The definition of the parsec is the distance which gives a parallax of one 
arcsecond, and is equivalent to 3.26 light-years, or 3.09 x 10 16 m. The field 
of parallax measurement was revolutionised by the Hipparcos satellite, which 
measured thousands of stellar parallax distances, including observations of 223 
Galactic Cepheids; of the Cepheids, 26 yielded determinations of reasonable 
significance [ITj . 

Some relatively nearby stars exist in clusters of a few hundred stars known as 
"open clusters" . These stars can be plotted on a Hertzsprung-Russell diagram 
of temperature, deduced from their colour together with Wien's law, against 
apparent luminosity. Such plots reveal a characteristic sequence, known as the 
"main sequence" which ranges from red, faint stars to blue, bright stars. This 
sequence corresponds to the main phase of stellar evolution which stars occupy 
for most of their lives when they are stably burning hydrogen. In some nearby 
clusters, notably the Hyades, we have stars all at the same distance and for 
which parallax effects can give the absolute distance. In such cases, the main 
sequence can be calibrated so that we can predict the absolute luminosity of a 
main-sequence star of a given colour. Applying this to other clusters, a process 
known as "main sequence fitting" , can also give the absolute distance to these 
other clusters. 

The next stage of the bootstrap process is to determine the distance to 
the nearest objects outside our own Galaxy, the Large and Small Magellanic 
Clouds. For this we can apply the open-cluster method directly, by observing 
open clusters in the LMC. Alternatively, we can use calibrators whose true 
luminosity we know, or can predict from their other properties. Such calibrators 
must be present in the LMC and also in open clusters (or must be close enough 
for their parallaxes to be directly measurable). 

These calibrators include Mira variables, RR Lyrae stars and Cepheid vari- 
able stars, of which Cepheids arc intrinsically the most luminous. All of these 
have variability periods which are correlated with their absolute luminosity, and 
in principle the measurement of the distance of a nearby object of any of these 
types can then be used to determine distances to more distant similar objects 
simply by observing and comparing the variability periods. 

The LMC lies at about 50 kpc, about three orders of magnitude less than 
that of the distant galaxies of interest for the Hubble constant. However, one 
class of variable stars, Cepheid variables, can be seen in both the LMC and 
in galaxies at much greater distances. The coming of the Hubble Space Tele- 
scope has been vital for this process, as only with the HST can Cepheids be 
reliably identified and measured in such galaxies. It is impossible to overstate 
the importance of Cepheids; without them, the connection between the LMC 
and external galaxies is very hard to make. 

Even the HST galaxies containing Cepheids are not sufficient to allow the 
measurement of the universal expansion, because they are not distant enough 
for the dominant velocity to be the Hubble flow. The final stage is to use 
galaxies with distances measured with Cepheid variables to calibrate other indi- 
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cators which can be measured to cosmologically interesting distances. The most 
promising indicator consists of type la supernovae (SNe), which are produced 
by binary systems in which a giant star is dumping mass on to a white dwarf 
which has already gone through its evolutionary process and collapsed to an 
electron-degenerate remnant; at a critical point, the rate and amount of mass 
dumping is sufficient to trigger a supernova explosion. The physics of the explo- 
sion, and hence the observed light-curve of the rise and slow fall, has the same 
characteristic regardless of distance. Although the absolute luminosity of the ex- 
plosion is not constant, type la supernovae have similar light-curves |119l l51 [T15] 
and in particular there is a very good correlation between the peak brightness 
and the degree of fading of the supernova 15 day:fl after peak brightness (a 
quantity known as Ami5 [1181 [53] ). If SNe la can be detected in galaxies with 
known Cepheid distances, this correlation can be calibrated and used to de- 
termine distances to any other galaxy in which a SN la is detected. Because 
of the brightness of supernovae, they can be observed at large distances and 
hence, finally, a comparison between redshift and distance will give a value of 
the Hubble constant. 

There are alternative indicators which can be used instead of SNe la for de- 
termination of H ; all of them rely on the correlation of some easily observable 
property of galaxies with their luminosity. For example, the edge-on rotation 
velocity v of spiral galaxies scales with luminosity as L cx v 4 , a scaling known 
as the Tully-Fisher relation [162J. There is an equivalent for elliptical galaxies, 
known as the Faber-Jackson relation [36j . In practice, more complex combi- 
nations of observed properties are often used such as the D n parameter of [34] 
and [91] , to generate measurable properties of elliptical galaxies which correlate 
well with luminosity, or the "fundamental plane" [34] [32] between three prop- 
erties, the luminosity within an effective radiu^l, the effective radius, and the 
central stellar velocity dispersion. Here again, the last two parameters are mea- 
surable. Finally, the degree to which stars within galaxies are resolved depends 
on distance, in the sense that closer galaxies have more statistical "bumpiness" 
in the surface-brightness distribution [157] , This method of surface brightness 
fluctuation can again be calibrated by Cepheid variables in the nearer galaxies. 

2.3 Problems and comments 
2.3.1 Distance to the LMC 

The LMC distance is probably the best-known, and least controversial, part of 
the distance ladder. Some methods of determination are summarised in [40] 
and little has changed since then. Independent calibrations using RR Lyrae 
variables, Cepheids and open clusters, are consistent with a distance of ~ 50 kpc. 
While all individual methods have possible systematics (see in particular the 

4 Because of the expansion of the Universe, there is a time dilation of a factor (1 + z)~ x 
which must be applied to timescales measured at cosmological distances before these are used 
for such comparisons. 

5 Thc effective radius is the radius from within which half the galaxy's light is emitted. 
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next Section T2.3.2I in the case of Cepheids), their agreement within the errors 
leaves little doubt that the measurement is correct. Moreover, an independent 
measurement was made in [108] using the type II supernova 1987A in the LMC. 
This supernova produced an expanding ring whose angular diameter could be 
measured using the HST. An absolute size for the ring could also be deduced 
by monitoring ultraviolet emission lines in the ring and using light travel time 
arguments, and the distance of 51.2 ± 3.1 kpc followed from comparison of the 
two. An extension to this light-echo method was proposed in |144j which exploits 
the fact that the maximum in polarization in scattered light is obtained when 
the scattering angle is 90°. Hence, if a supernova light echo were observed in 
polarized light, its distance would be unambiguously calculated by comparing 
the light-echo time and the angular radius of the polarized ring. 

The distance to the LMC adopted by most researchers in the field is between 
fi° = 18.50 and 18.54, in the units of "distance modulus" (defined as 51ogd— 5, 
where d is the distance in parsecs) corresponding to a distance of 50-51 kpc. 
The likely error in Ho of ~ 2% is well below the level of systematic errors in 
other parts of the distance ladder; recent developments in the use of standard- 
candle stars, main sequence fitting and the details of SN 1987A are reviewed 
in [5] where it is concluded that ^° = 18.50 ± 0.02. 

2.3.2 Cepheids 

If the Cepheid period-luminosity relation were perfectly linear and perfectly 
universal (that is, if we could be sure that it applied in all galaxies and all 
environments) the problem of transferring the LMC distance outwards to exter- 
nal galaxies would be simple. Unfortunately, to very high accuracy it may be 
neither. Although there are other systematic difficulties in the distance ladder 
determinations, problems involving the physics and phenomenology of Cepheids 
are currently the most controversial part of the error budget, and are the pri- 
mary source of differences in the derived values of H . 

The largest samples of Cepheids outside our own Galaxy come from mi- 
crolensing surveys of the LMC, reported in |164j . Sandage et al. [133j reanalyse 
those data for LMC Cepheids and claim that the best fit involves a break in the 
P-L relation at P ~ 10 days. In all three HST colours (B, V, I) the resulting 
slopes are different from the Galactic slopes, in the sense that at long periods, 
Galactic Cepheids are brighter than LMC Cepheids and are fainter at short pe- 
riods. The period at which LMC and Galactic Cepheids have equal luminosities 
is approximately 30 days in B, but is a little more than 10 days in i0. Sandage 
et al. [133J therefore claim a colour-dependent difference in the P-L relation 
which points to an underlying physical explanation. The problem is potentially 
serious in that the difference between Galactic and LMC Cepheid brightness can 
reach 0.3 magnitudes, corresponding to a 15% difference in inferred distance. 

6 Nearly all Cepheids measured in galaxies containing a SN la have periods > 20 days, so 
the usual sense of the effect is that Galactic Cepheids of a given period are brighter than LMC 
Cepheids. 
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At least part of this difference is almost certainly due to metallicity ef- 
fect!]. Groenewegen et al. [51] assemble earlier spectroscopic estimates of 
metallicity in Cepheids both from the Galaxy and the LMC and compare them 
with their independently derived distances, obtaining a marginally significant 
(—0.8 ± 0.3 mag dex -1 ) correlation of brightness with increasing metallicity by 
using only Galactic Cepheids. Using also the LMC cepheids gives —0.27 ± 
0.08 mag dex -1 . 

In some cases, independent distances to galaxies are available in the form 
of studies of the tip of the red giant branch. This phenomenon refers to the 
fact that metal-poor, population II red giant stars have a well-defined cutoff in 
luminosity which, in the /-band, does not vary much with nuisance parameters 
such as stellar population age. Deep imaging can therefore provide an indepen- 
dent standard candle which can be compared with that of the Cepheids, and 
in particular with the metallicity of the Cepheids in different galaxies. The re- 
sult [130] is again that metal-rich Cepheids are brighter, with a quoted slope of 
—0.24 ± 0.05 mag dex -1 . This agrees with earlier determinations [75] IZ2 and 
is usually adopted when a global correction is applied. 

The LMC is relatively metal-poor compared to the Galaxy, and the same 
appears to be true of its Cepheids. On average, the Galactic Cepheids tabulated 
in [ST] are approximately of solar metallicity, whereas those of the LMC are 
approximately —0.6 dex less metallic, corresponding to an 8% distance error 
if no correction is applied in the bootstrapping of Galactic to LMC distance. 
Hence, a metallicity correction must be applied when using the highest quality 
P-L relations from the OGLE observations of LMC Cepheids to the typically 
more metallic Cepheids in galaxies with SNe la observations. 

2.4 Independent local distance-scale methods 
2.4.1 Masers 

One exception to the rule that Cepheids are necessary for tying local and more 
global distance determinations is provided by the study of masers, the prototype 
of which is the maser system in the galaxy NGC 4258 [53]. This galaxy has 
a shell of masers which are oriented almost edge-on [961 150] and apparently in 
Keplerian rotation. As well as allowing a measurement of the mass of the central 
black hole, the velocity drift (acceleration) of the maser lines from individual 
maser features can also be measured. This allows a measurement of absolute 
size of the maser shell, and hence the distance to the galaxy. This has become 
steadily more accurate since the original work [54] [66] [3] . Macri et al. [92] also 
measure Cepheids in this object to determine a Cepheid distance (see Figure [3]) 
and obtain consistency with the maser distance provided that the LMC distance, 
to which the Cepheid scale is calibrated, is 48 ± 2 kpc. 

Further discoveries and observations of masers could in principle establish a 

7 Here, as elsewhere in astronomy, the term "metals" is used to refer to any element heavier 
than helium. Metallicity is usually quoted as 12+log(0 /H), where O and H are the abundances 
of oxygen and hydrogen. 



12 



S40S5 p. 9.16 strsa p. is.ae 




0.5 1 1.5 2 0.5 1 1.5 



Figure 3: Positions of Cepheid variables in HST/ACS observations of the galaxy 
NGC 4258, reproduced from [92] (upper panel) . Typical Cepheid lightcurves are 
shown in the lower panel. 
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distance ladder without heavy reliance on Cepheids. The Water Maser Cosmol- 
ogy Project 

( jhttp : //www . cf a . harvard . edu/wmcp/ index . html[ ) is conducting monitoring 
and high-resolution imaging of samples of extragalactic masers, with the even- 
tual aim of a maser distance scale accurate to ~ 3%. 

2.4.2 Other methods of establishing the distance scale 

Several other different methods have been proposed to bypass some of the early 
rungs of the distance scale and provide direct measurements of distance to rela- 
tively nearby galaxies. Many of these are reviewed in the article by Oiling [105] . 

One of the most promising methods is the use of detached eclipsing binary 
stars to determine distances directly [107 . In nearby binary stars, where the 
components can be resolved, the determination of the angular separation, period 
and radial velocity amplitude immediately yields a distance estimate. In more 
distant eclipsing binaries in other galaxies, the angular separation cannot be 
measured directly. However, the light-curve shapes provide information about 
the orbital period, the ratio of the radius of each star to the orbital separation, 
and the ratio of the stars' luminosities. Radial velocity curves can then be used 
to derive the stellar radii directly. If we can obtain a physical handle on the 
stellar surface brightness (e.g. by study of the spectral lines) then this, together 
with knowledge of the stellar radius and of the observed flux received from each 
star, gives a determination of distance. The DIRECT project [H] has used 
this method to derive a distance of 964 ± 54 kpc to M33, which is higher than 
standard distances of 800-850 kpc :46, 87]. It will be interesting to see whether 
this discrepancy continues after further investigation. 

A somewhat related method, but involving rotations of stars around the 
centre of a distant galaxy, is the method of rotational parallax [1171 11061 1105] . 
Here one observes both the proper motion corresponding to circular rotation, 
and the radial velocity, of stars within the galaxy. Accurate measurement of 
the proper motion is difficult and will require observations from future space 
missions. 

2.5 H : 62 or 73? 

We are now ready to try to disentangle and understand the reasons why inde- 
pendent analyses of the same data give values which are discrepant by twice the 
quoted systematic errors. Probably the fairest and most up-to-date analysis is 
achieved by comparing the result of Riess et al. |124 i] from 2005 (R05) who found 
Hq = 73 km s _1 Mpc -1 , with statistical errors of 4 km s _1 Mpc -1 , and the one 
of Sandage et al. [TMj from 2006 (S06) who found 62.3±1.3 km s -1 Mpc -1 (sta- 
tistical). Both papers quote a systematic error of 5 km s -1 Mpc -1 . The R05 
analysis is based on four SNe la: 1994ae in NGC 3370, 1998aq in NGC 3982, 
1990N in NGC 4639 and 1981B in NGC 4536. The S06 analysis includes 
eight other calibrators, but this is not an issue as S06 find Ho = 63.3 ± 
1.9 km s -1 Mpc -1 from these four calibrators separately, still a 15% difference. 
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Inspection of Table 13 of R05 and Table 1 of S06 reveals part of the problem; 
the distances of the four calibrators are generally discrepant in the two analyses, 
with S06 having the higher value. In the best case, SN 1994ae in NGC 3370, the 
discrepancy is only 0.08 ma£|. In the worst case (SN 1990N in NGC 4639) R05 
quotes a distance modulus fjP = 31.74, whereas the value obtained by S06 is 
32.20, corresponding to a 20-25% difference in the inferred distance and hence 
in Hq. The quoted fir is formed by a combination of observations at two optical 
bands, V and /, and is normally defined as 2.52(m — M)j — 1.52(m — M)y, 
although the coefficients differ slightly between different authors. The purpose 
of the combination is to eliminate differential effects due to reddening, which 
it does exactly provided that the reddening law is known. This law has one 
parameter, R, known as the "ratio of total to selective extinction", and defined 
as the number of magnitudes of extinction at V corresponding to one magnitude 
of difference between B and V. 

We can investigate what is going on if we go back to the original photometry. 
This is given in [126] and has been corrected for various effects in the WFPC2 
photometry discovered since the original work of Holtzman et al. 57jf|. If we 
follow R05, we proceed by using the Cepheid P-L relation for the LMC given 
by [156] . We then apply a global metallicity correction, to account for the fact 
that the LMC is less metallic than NGC 4639 by about 0.6 dex [TSD], and we 
arrive at the R05 value for fi°. Alternatively, we can use the P-L relations given 
by S06 and derived from earlier work by Sandage et al. [133] . These authors 
derived relations both for the LMC and for the Galaxy. Like the P-L relations 
in [156] . the LMC relations are based on the OGLE observations in [164] . with 
the addition of further Cepheids at long periods, and the Galactic relations are 
based on earlier work on Galactic cepheids EU EH [73]. Following S06, 
we derive /x° values separately for each Cepheid for the Galaxy and LMC. We 
then assume that the difference between the LMC and Galactic P-L relations 
is entirely due to metallicity, and use the measured NGC 4639 metallicity to 
interpolate and find the corrected n° value for each Cepheid. This interpolation 
gives us much larger metallicity corrections than in R05. We then finally average 
the fjP values to recover S06's longer distance modulus. 

The major difference is not specifically in the P-L relation assumed for the 
LMC, because the relation in [156j used by R05 is virtually identical to the 
P-L relation for long-period Cepheids used by S06. The difference lies in the 
correction for metallicity. R05 use a global correction 

Amr 05 = 0.24A[O/H] (10) 

from [130] (A[0/H] is the metallicity of the observed Cepheids minus the metal- 
licity of the LMC), whereas S06's correction by interpolation between LMC and 

8 Indeed, R05 calculate the value of Ho for SN 1994ae together with SN 1998aq according 
to the prescription of the Sandage et al. group as at 2005, and find 69 km s" 1 Mpc - 1 . 

9 There are two effects here. The first is the "long versus short" effect, which causes a 
decrease of recorded flux of a few percent in short exposures compared to long ones. The 
second is the effect of radiation damage, which affected later WFPC2 observations more than 
earlier ones and resulted in a uniform decrease of charge transfer efficiency and observed flux. 
This is again an effect at the few percent level. 
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Galactic P-L relations is [126] 



A MS06 = 1-67 (logP - 0.933)A[O/H]. (11) 

Which is correct? Both methods are perfectly defensible given the assump- 
tions that are made. The S06 crucially depends on the Galactic P-L relations 
being correct, and in addition depends on the hypothesis that the difference 
between Galactic and LMC P-L relations is dominated by metallicity effects 
(although it is actually quite hard to think of other effects that could have any- 
thing like the same systematic effect as the composition of the stars involved). 
The S06 Galactic relations are based on Tammann et al. [150 (TSR03) who 
in turn derive them from two other methods. The first is the calibration of 
Cepheids in open clusters to which the distance can be measured independently 
(see Section [2~ljl . as applied in [40]. The second is a compilation in [48] including 
earlier measurements and compilations (see e.g. [43]) of stellar angular diame- 
ters by lunar occultation and other methods. Knowing the angular diameters 
and temperatures of the stars, distances can be determined [1711 [5] essentially 
from Stefan's law. These two methods are found to agree in |150j . but this 
agreement and the consequent steep P-L relations for Galactic Cepheids, are 
crucial to the S06 case. Macri et al. [92] explicitly consider this assumption 
using new ACS observations of Cepheids in the galaxy NGC 4258 which have 
a range of metallicities [179] . They find that, if they assume the P-L relations 
of TSR03 whose slope varies with metallicity, the resulting /jP determined from 
each Cepheid individually varies with period, suggesting that the TSR03 P-L re- 
lation overcorrects at long period and hence that the P-L assumptions of R05 
are more plausible. It is probably fair to say that more data is needed in this 
area before final judgements are made. 

The R05 method relies only on the very well-determined OGLE photometry 
of the LMC Cepheids and not on the Galactic measurements, but does rely on a 
global metallicity correction which is vulnerable to period-dependent metallicity 
effects. This is especially true since the bright Cepheids typically observed in 
relatively distant galaxies which host SNe la are weighted towards long periods, 
for which S06 claim that the metallicity correction is much larger. 

Although the period-dependent metallicity correction is a major effect, there 
are a number of other differences which affect Ho by a few percent each. 

The calibration of the type la supernova distance scale, and hence H , is 
affected by the selection of galaxies used which contain both Cepheids and 
historical supernovae. Riess et al. [124] make the case for the exclusion of a 
number of older supernovae with measurements on photographic plates. Their 
exclusion, leaving four calibrators with data judged to be of high quality, has the 
effect of shrinking the average distances, and hence raising Hq, by a few percent. 
Freedman et al. [45] included six galaxies including SN 1937C, excluded in [124] . 
but obtained approximately the same value for Hq. 

There is a selection bias in Cepheid variable studies in that faint Cepheids are 
harder to see. Combined with the correlation between luminosity and period, 
this means that only the brighter short-period Cepheids are seen, and therefore 
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that the P-L relation in distant galaxies is made artificially shallow [132] result- 
ing in underestimates of distances. Neglect of this bias can give differences of 
several percent in the answer, and detailed simulations of it have been carried 
out by Teerikorpi and collaborators (see e.g. |152Ulllllll2U113] ). Most authors 
correct explicitly for this problem - for example, Freedman et al. [15] calculate 
the correction analytically and find a maximum bias of about 3%. Teerikorpi 
and Paturel suggest that a residual bias may still be present, essentially because 
the amplitude of variation introduces an additional scatter in brightness at a 
given period, in addition to the scatter in intrinsic luminosity. How big this bias 
is is hard to quantify, although it can in principle be eliminated by using only 
long-period Cepheids at the cost of increases in the random error. 

Further possible effects include differences in SNe la luminosities as a func- 
tion of environment. Wang et al. [170) used a sample of 109 supernovae to 
determine a possible effect of metallicity on SNe la luminosity, in the sense that 
supernovae closer to the centre of the galaxy (and hence of higher metallic- 
ity) are brighter. They include colour information using the indicator AC12 = 
(B — V)i2 days, the B — V colour at 12 days after maximum, as a means of re- 
ducing scatter in the relation between peak luminosity and Amis which forms 
the traditional standard candle. Their value of H is, however, quite close to the 
Key Project value, as they use the four galaxies of [124] to tie the supernova and 
Cepheid scales together. This closeness indicates that the SNe la environment 
dependence is probably a small effect compared with the systematics associated 
with Cepheid metallicity 

In summary, local distance measures have converged to within 15%, a vast 
improvement on the factor of 2 uncertainty which prevailed until the late 1980s. 
Further improvements are possible, but involve the understanding of some non- 
trivial systematics and in particular require general agreement on the physics 
of metallicity effects on Cepheid P-L relations. 
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3 The CMB and Cosmological Estimates of the 
Distance Scale 



3.1 The physics of the anisotropy spectrum and its impli- 
cations 

The physics of stellar distance calibrators is very complicated, because it comes 
from the era in which the Universe has had time to evolve complicated astro- 
physics. A large class of alternative approaches to cosmological parameters in 
general involve going back to a substantially astrophysics-free zone, the epoch 
of recombination. Although none of these tests uniquely or directly determine 
Hq, they provide joint information about Hq and other cosmological parameters 
which is improving at a very rapid rate. 

In the Universe's early history, its temperature was high enough to prohibit 
the formation of atoms, and the Universe was therefore ionized. Approximately 
3 x 10 5 yr after the Big Bang, corresponding to a redshift z rec ~ 1000, the 
temperature dropped enough to allow the formation of atoms, a point known 
as "recombination" . For photons, the consequence of recombination was that 
photons no longer scattered from ionized particles but were free to stream. After 
recombination, these primordial photons reddened with the expansion of the 
Universe, forming the cosmic microwave background (CMB) which we observe 
today as a black-body radiation background at 2.73 K. 

In the early Universe, structure existed in the form of small density fluctua- 
tions (Sp/ p ~ 0.01) in the photon-baryon fluid. The resulting pressure gradients, 
together with gravitational restoring forces, drove oscillations, very similar to 
the acoustic oscillations commonly known as sound waves. At the same time, 
the Universe expanded until recombination. At this point, the structure was 
dominated by those oscillation frequencies which had completed a half-integral 
number of oscillations within the characteristic size of the Universe at recombi- 
nation; this pattern became frozen into the photon field which formed the CMB 
once the photons and baryons decoupled. The process is reviewed in much more 
detail in [50] , 

The resulting "acoustic peaks" dominate the fluctuation spectrum (see Fig- 
ure |4|) . Their angular scale is a function of the size of the Universe at the time 
of recombination, and the angular diameter distance between us and z rcc . Since 
the angular diameter distance is a function of cosmological parameters, measure- 
ment of the positions of the acoustic peaks provides a constraint on cosmological 
parameters. Specifically, the more closed the spatial geometry of the Universe, 
the smaller the angular diameter distance for a given redshift, and the larger 
the characteristic scale of the acoustic peaks. The measurement of the peak 
position has become a strong constraint in successive observations (in particu- 
lar Boomerang, reported in [3D] and WMAP, reported in [146] and |145j ) and 
corresponds to an approximately spatially flat Universe in which Q m + fl\ ~ 1. 

But the global geometry of the Universe is not the only property which 
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Figure 4: Diagram of the CMB anisotropics, plotted as strength against spatial 
frequency, from the WMAP 3-year data [1451. The measured points are shown 
together with best-fit models to the 1-year and 3-year WMAP data. Note the 
acoustic peaks, the largest of which corresponds to an angular scale of about half 
a degree. 

can be deduced from the fluctuation spectrum^. The peaks are also sensitive 
to the density of baryons, of total (baryonic plus dark) matter, and of dark 
energy (energy associated with the cosmological constant or more generally 
with w < — | components) . These densities scale with the square of the Hubble 
parameter times the corresponding dimensionless densities (see Equation |5jl) 
and measurement of the acoustic peaks therefore provides information on the 
Hubble constant, degenerate with other parameters, principally the curvature 
energy ilk and the index w in the dark energy equation of state. The second peak 
strongly constrains the baryon density, H^Hq, and the third peak is sensitive to 
the total matter density in the form Cl m Hg. 

3.2 Degeneracies and implications for H 

If the Universe is exactly flat and the dark energy is a cosmological constant, 
then the debate about Hq is over. The Wilkinson Anisotropy Probe (WMAP) 
has so far published two sets of cosmological measurements, based on one year 
and three years of operation [1461 1145j . From the latest data, the assumption 
of a spatially flat Universe requires H = 73 ± 3 km s _1 Mpc , and also deter- 
mines other cosmological parameters to two and in some cases three significant 
figures [145] . If we do not assume the Universe to be exactly flat, then we obtain 
a degeneracy with Hq in the sense that every decrease of 20 km s _1 Mpc -1 in- 
creases the total density of the Universe by 0.1 in units of the closure density (see 

10 See http : //background .uchicago . edu^whu/ intermediate/intermediate .html for a 

much longer exposition and tutorial on all these areas. 
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Figure [5]). The WMAP data by themselves, without any further assumptions 
or extra data, do not supply a significant constraint on Hq. 

There are two other major programmes which result in constraints on com- 
binations of H , flm, (now considered as a general density in dark energy 
rather than specifically a cosmological constant energy density) and w. The first 
is the study of type la supernovae, which as we have seen function as standard 
candles, or at least easily calibratable candles. Studies of supernovae at cosmo- 
logical redshifts by two different collaborations [1161 11151 1123] have shown that 
distant supernovae are fainter than expected if the simplest possible spatially 
flat model (the Einstein-de Sitter model, for which Q m = 1, Q\ = 0) is correct. 
The resulting determination of luminosity distance has given constraints in the 
^m _ ^A plane which are more or less orthogonal to the WMAP constraints. 

The second important programme is the measurement of structure at more 
recent epochs than the epoch of recombination. This is interesting because fluc- 
tuations prior to recombination can propagate at the relativistic (c/v3) sound 
speed which predominates at that time. After recombination, the sound speed 
drops, effectively freezing in a characteristic length scale to the structure of 
matter which corresponds to the propagation length of acoustic waves by the 
time of recombination. This is manifested in the real Universe by an expected 
preferred correlation length of ~ 100 Mpc between observed baryon structures, 
otherwise known as galaxies. The largest sample available for such studies comes 
from luminous red galaxies (LRGs) in the Sloan Digital Sky Survey [177] . The 
expected signal has been found |35j in the form of an increased power in the 
cross-correlation between galaxies at separations of about 100 Mpc. It corre- 
sponds to an effective measurement of angular diameter distance to a redshift 
z ~ 0.35. 

As well as supernova and acoustic oscillations, several other slightly less 
tightly-constraining measurements should be mentioned: 

• Lyman a forest observations. The spectra of distant quasars have deep 
absorption lines corresponding to absorbing matter along the line of sight. 
The distribution of these lines measures clustering of matter on small 
scales and thus carries cosmological information (see e.g. [163, 95J). 

• Clustering on small scales |153j . The matter power spectrum can be mea- 
sured using large samples of galaxies, giving constraints on combinations 
of Hq, Q m and erg, the normalization of matter fluctuations on 8-Mpc 
scales. 

3.2.1 Combined constraints 

Tegmark et al. |154j have considered the effect of applying the SDSS acoustic os- 
cillation detection together with WMAP data. As usual in these investigations, 
the tightness of the constraints depends on what is assumed about other cos- 
mological parameters. The maximum set of assumptions (called the "vanilla" 
model by Tegmark et al.) includes the assumption that the spatial geometry 
of the Universe is exactly flat, that the dark energy contribution results from 
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a pure w = — 1 cosmological constant, and that tensor modes and neutrinos 
make neglible contributions. Unsurprisingly, this gives them a very tight con- 
straint on Hq of 73 ± 1.9 km s _1 Mpc" 1 . However, even if we now allow the 
Universe to be not exactly flat, the use of the detection of baryon acoustic oscil- 
lations in [35) together with the WMAP data yields a 5%-error measurement of 
Hq = 7I.6I4 3 km s _1 Mpc -1 . This is entirely consistent with the Hubble Key 
Project measurement from Cepheid variables, but only just consistent with the 
version in [134] . The improvement comes from the extra distance measurement, 
which provides a second joint constraint on the variable set {Hq, f2 m , w). 

Even this value, however, makes the assumption that w — —1. If we relax 
this assumption as well, Tegmark ct al. [154] find that the constraints broaden 
considerably, to the point where the 2a bounds on Hq range lie between 61 
and 84 km s" 1 Mpc" 1 (see Figure [5] and [154] 1. even if the HST Key Project 
results [45] are added. It has to be said that both w = — 1 and Qk = are 
highly plausible assumption^]], and if only one of them is correct, Hq is known 
to high accuracy. To put it another way, however, an independent measurement 
of Hq would be extremely useful in constraining all of the other cosmological 
parameters provided that its errors were at the 5% level or better. In fact [55] . 
"The single most important complement to the CMB for measuring the dark 
energy equation of state at z > 0.5 is a determination of the Hubble constant to 
better than a few percent" . Oiling [105] quantifies this statement by modelling 
the effect of improved Hq estimates on the determination of w. He finds that, 
although a 10% error on Hq is not a significant contribution to the current 
error budget on w, but that once improved CMB measurements such as those 
to be provided by the Planck satellite are obtained, decreasing the errors on Hq 
by a factor of five to ten could have equal power to much of the (potentially 
more expensive, but in any case usefully confirmatory) direct measurements of 
w planned in the next decade. In the next Section [4] explore various ways by 
which this might be achieved. 

It is of course possible to put in extra information, at the cost of introducing 
more data sets and hence more potential for systematic error. Inclusion of the 
supernova data, as well as Ly-a forest data [p_5], SDSS and 2dF galaxy cluster- 
ing [HHHl] and other CMB experiments (CBI, [T2T3] : VSA, [31]; Boomerang, 
[93] . Acbar, [85]), together with a vanilla model, unsurprisingly gives a very 
tight constraint on Hq of 70.5 ± 1.3 km s _1 Mpc" 1 . Including non-vanilla pa- 
rameters one at a time also gives extremely tight constraints on the spatial 
flatness (O fe = -0.003 ± 0.006) and w (-1.04 ± 0.06), but the constraints are 
again likely to loosen if both w and £lk are allowed to depart from vanilla val- 
ues. A vast literature is quickly assembling on the consequences of shoehorning 
together all possible combinations of different datasets with different parameter 
assumptions (see e.g. [23 H ED [51 [T75J EU). 



xl Beware plausible assumptions, however; fifteen years ago Q\ = was a highly plausible 
assumption. 
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Figure 5: Top: The allowed range of the parameters f2 m , £Ia, from the WMAP 3- 
year data, is shown as a series of points (reproduced from 114-51). The diagonal 
line shows the locus corresponding to a flat Universe (T2 m + = 1), An 
exactly flat Universe corresponds to Hq ~ 70 km s - Mpc^ 1 , but lower values 
are allowed provided the Universe is slightly closed. Bottom: Analysis reproduced 
from 1 154*1 showing the allowed range of the Hubble constant, in the form of 
the Hubble parameter h = Hq/100 km s - Mpc - 1 , and Cl m by combination of 
WMAP 3-year data with acoustic oscillations. A range of Hq is still allowed by 
these data, although the allowed region shrinks considerably if we assume that 
w = — 1 or Slfc = 0. 
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4 One-Step Distance Methods 



4.1 Gravitational lenses 

A general review of gravitational lensing is given in [169]; here we review the 
theory necessary for an understanding of the use of lenses in determining the 
Hubble constant. 

4.1.1 Basics of lensing 

Light is bent by the action of a gravitational field. In the case where a galaxy 
lies close to the line of sight to a background quasar, the quasar's light may 
travel along several different paths to the observer, resulting in more than one 
image. 

The easiest way to visualise this is to begin with a zero-mass galaxy (which 
bends no light rays) acting as the lens, and considering all possible light paths 
from the quasar to the observer which have a bend in the lens plane. From 
the observer's point of view, we can connect all paths which take the same 
time to reach the observer with a contour, which in this case is circular in 
shape. The image will form at the centre of the diagram, surrounded by circles 
representing increasing light travel times. This is of course an application of 
Fermat's principle; images form at stationary points in the Fermat surface, in 
this case at the Fermat minimum. Put less technically, the light has taken a 
straight-line path between the source and observer. 

If we now allow the galaxy to have a steadily increasing mass, we introduce 
an extra time delay (known as the Shapiro delay) along light paths which pass 
through the lens plane close to the galaxy centre. This makes a distortion in the 
Fermat surface. At first, its only effect is to displace the Fermat minimum away 
from the distortion. Eventually, however, the distortion becomes big enough to 
produce a maximum at the position of the galaxy, together with a saddle point 
on the other side of the galaxy from the minimum. By Fermat's principle, two 
further images will appear at these two stationary points in the Fermat surface. 
This is the basic three-image lens configuration, although in practice the central 
image at the Fermat maximum is highly demagnified and not usually seen. 

If the lens is significantly elliptical and the lines of sight are well aligned, 
we can produce five images, consisting of four images around a ring alternating 
between maxima and saddle points, and a central, highly demagnified Fermat 
maximum. Both four- image and two- image systems ("quads" and "doubles") 
are in fact seen in practice. The major use of lens systems is for determining 
mass distributions in the lens galaxy, since the positions and brightnesses of the 
images carry information about the gravitational potential of the lens. Gravi- 
tational lensing has the advantage that its effects are independent of whether 
the matter is light or dark, so in principle the effects of both baryonic and 
non-baryonic matter can be probed. 
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4.1.2 Principles of time delays 



Refsdal [122] pointed out that if the background source is variable, it is possible 
to measure an absolute distance within the system and therefore the Hubble 
constant. To see how this works, consider the light paths from the source to 
the observer corresponding to the individual lensed images. Although each is 
at a stationary point in the Fermat time delay surface, the absolute light travel 
time for each will generally be different, with one of the Fermat minima having 
the smallest travel time. Therefore, if the source brightens, this brightening will 
reach the observer at different times corresponding to the two different light 
paths. Measurement of the time delay corresponds to measuring the difference 
in the light travel times, each of which is individually given by 



where a, (3 and 9 are angles defined above in Figure D\, D s and D\ s are 
angular diameter distances also defined in Figure [H z\ is the lens redshift, and 
ip(6) is a term representing the Shapiro delay of light passing through a gravi- 
tational field. Fermat's principle corresponds to the requirement that Vr = 0. 
Once the differential time delays are known, we can then calculate the ratio 
of angular diameter distances which appears in the above equation. If the 
source and lens redshifts are known, Hq follows immediately. A handy rule of 
thumb which can be derived from this equation for the case of a 2-image lens, 
if we make the assumption that the matter distribution is isothermal 12 ! and 
H a = 70 km s _1 Mpc -1 , is 



where z\ is the lens redshift, s is the separation of the images (approximately 
twice the Einstein radius) , / > 1 is the ratio of the fluxes and D is the value of 
D S D\/ D\ s in Gpc. A larger time delay implies a correspondingly lower Hq. 

The first gravitational lens was discovered in 1979 [168] and monitoring pro- 
grammes began soon afterwards to determine the time delay. This turned out to 
be a long process involving a dispute between proponents of a ~ 400— day and 
a ~ 550— day delay, and ended with a determination of 417 ± 2 days [541, 1136] . 
Since that time, 17 more time delays have been determined (see Table [1]). In 
the early days, many of the time delays were measured at radio wavelengths 
by examination of those systems in which a radio-loud quasar was the multiply 
imaged source (see Figure [7]) . Recently, optically- measured delays have domi- 
nated, due to the fact that only a small optical telescope in a site with good 
seeing is needed for the photometric monitoring, whereas radio time delays re- 
quire large amounts of time on long-baseline interferometers which do not exist 
in large numberj^l. 

12 An isothermal model is one in which the projected surface mass density decreases as 1/r. 
An isothermal galaxy will have a flat rotation curve, as is observed in many galaxies. 

13 Essentially all radio time delays have come from the VLA, although monitoring pro- 
grammes with MERLIN have also been attempted. 




(12) 




(13) 
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Figure 6: Basic geometry of a gravitational lens system, reproduced from [169] . 
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Figure 7: The lens system JVAS B0218+357. On the right is shown the mea- 
surement of time delay of about 10 days from asynchronous variations of the 
two lensed images [10]. The left panels show the HST/ACS image \178^ on 
which can be seen the two images and the spiral lensing galaxy, and the radio 
MERLIN+VLA image Jllf showing the two images together with an Einstein 
ring. 

4.1.3 The problem with lens time delays 

Unlike local distance determinations (and even unlike cosmological probes which 
typically use more than one measurement), there is only one major systematic 
piece of astrophysics in the determination of H by lenses, but it is a very 
important one. This is the form of the potential in Equation ([T2"|) . If one 
parametrises the potential in the form of a power law in projected mass density 
versus radius, the index is —1 for an isothermal model. This index has a pretty 
direct degeneracj[^3 with the deduced length scale and therefore the Hubble 
constant; for a change of 0.1, the length scale changes by about 10%. The 
sense of the effect is that a steeper index, which corresponds to a more centrally 
concentrated mass distribution, decreases all the length scales and therefore 
implies a higher Hubble constant for a given time delay. The index cannot be 
varied at will, given that galaxies consist of dark matter potential wells into 
which baryons have collapsed and formed stars. The basic physics means that 
it is almost certain that matter cannot be more centrally condensed than the 
stars, and cannot be less centrally condensed than the theoretically favoured 
"universal" dark matter profile, known as a NFW profile jlOOj . 

14 As discussed extensively in | 77l I80 |, this is not a global degeneracy, but arises because 
the lensed images tell you about the mass distribution in the annulus centred on the galaxy 
and with inner and outer radii defined by the inner and outer images. Kochanek |77| derives 
detailed expressions for the time delay in terms of the central underlying and controlling 
parameter, the surface density in this annulus [49] . 
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Worse still, all matter along the line of sight contributes to the lensing po- 
tential in a particularly unpleasant way; if one has a uniform mass sheet in the 
region, it does not affect the image positions and fluxes which form the con- 
straints on the lensing potential, but it does affect the time delay. It operates in 
the sense that, if a mass sheet is present which is not known about, the length 
scale obtained is too short and consequently the derived value of H is too high. 
This mass-sheet degeneracy [49] can only be broken by lensing observations 
alone for a lens system which has sources at multiple redshifts, since there are 
then multiple measurements of angular diameter distance which are only con- 
sistent, for a given mass sheet, with a single value of Hq. Existing galaxy lenses 
do not contain examples of this phenomenon. 

Even worse still, there is no guarantee that parametric models describe lens 
mass distributions to the required accuracy. In a series of papers [1281 11721 11291 
1127] non-parametric, pixellated models of galaxy mass distributions have been 
developed which do not require any parametric form, but only basic physical 
plausibility arguments such as monotonic outwards decrease of mass density. 
Not surprisingly, error bars obtained by this method are larger than for para- 
metric methods, usually by factors of at least 2. 

4.1.4 Now and onwards in time delays and modelling 

Table Q] shows the currently measured time delays, with references and com- 
ments. Since the most recent review |80) an extra half-dozen have been added, 
and there is every reason to suppose that the sample will continue to grow at a 
similar rateF^I. 

Despite the apparently depressing picture painted in the previous Section r4.1.3l 
about the prospects for obtaining mass models from lenses, the measurement of 
Hq is improving in a number of ways. 

First, some lenses have more constraints on the mass model than others. 
The word "more" here is somewhat subjective, but examples include JVAS 
B0218+357 which in addition to two images, also has VLBI substructure within 
each image and an Einstein ring formed from an extended background source, 
and CLASS B1933+503 which has three background radio sources, each multi- 
ply imaged. Something of a Murphy's Law operates in the latter case, however, 
as compact triple radio sources tend to be of the class known as Compact Sym- 
metric Objects (CSOs) which do not vary and therefore do not give time delays. 
Einstein rings in general give good constraints [78j although non-parametric 
models are capable of soaking up most of these in extra degrees of freedom [129] . 
In general however, no "golden" lenses with multiple constraints and no other 
modelling problems have been founcP^l. The best models of all come from lenses 

15 The aspiring measurer of time delays is faced with a dilemma, in terms of whether to 
justify the proposal in terms of measuring Ho, given the previously mentioned problems with 
mass modelling, or in terms of determining mass models by assuming Hq = 71 km s _1 Mpc -1 
(or whatever). 

le This term, "golden lens", was coined over ten years ago to describe the perfect time-delay 
lens. Since none have been found, it is probably time to abandon it and metallic derivatives 
thereof, perhaps in return for vanilla CMB models and their associated error ellipses, vanilla 
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Table 1: Time delays, with 1-a errors, from the literature. In some cases mul- 
tiple delays have been measured in 4-image lens systems, and in this case each 
delay is given separately for the two components in brackets. 
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from the SLACS survey, which have extended sources [H] but unfortunately the 
previous Murphy's law applies here too; extended sources are not variable. 

Second, it is possible to increase the reliability of individual lens mass models 
by gathering extra information. A major improvement is available by the use 
of stellar velocity dispersions [1591 11581 11601 152"] measured in the lensing galaxy. 
As a standalone determinant of mass models in galaxies at z ~ 0.5, typical of 
lens galaxies, such measurements are not very useful as they suffer from severe 
degeneracies with the structure of stellar orbits. However, the combination of 
lensing information (which gives a very accurate measurement of mass enclosed 
by the Einstein radius) and stellar dynamics (which gives, more or less, the mass 
enclosed within the effective radius of the stellar light) gives a measurement that 
is in principle a very direct constraint on the mass slope. The method has large 
error bars, in part due to residual dependencies on the shape of stellar orbits, 
but also because these measurements are very difficult; each galaxy requires 
about one night of good seeing on a 10-m telescope. It is also not certain that 
the mass slope between Einstein and effective radii is always a good indicator 
of the mass slope within the annulus between the lensed images. Nevertheless, 
this programme has the extremely beneficial effect of turning a systematic error 
in each galaxy into a smaller, more-or-less random error. 

Third, we can remove problems associated with mass sheets associated with 
nearby groups by measuring them using detailed studies of the environments of 
lens galaxies. Recent studies of lens groups [55J ES GUI [57] show that neglecting 
matter along the line of sight typically has an effect of 10-20%, with matter 
close to the redshift of the lens contributing most. 

Finally, we can improve measurements in individual lenses which have mea- 
surement difficulties. For example, in the lenses 1830-211 [90] and B0218+357 [P09] 
the lens galaxy position is not well known. In the case of B0218+357, York 
et al. [178 present deep HST/ACS data which allow much better astrometry. 
Overall, by a lot of hard work using all methods together, the systematic errors 
involved in the mass model in each lens individually can be reduced to a random 
error. We can then study lenses as an ensemble. 

Early indications, using systems with isolated lens galaxies in uncomplicated 
environments, and fitting isothermal mass profiles, resulted in rather low values 
of the Hubble constant (in the high fifties [76]). In order to be consistent 
with Ho ~ 70 km s _1 Mpc -1 , the mass profiles had to be steepened to the 
point where mass followed light; although not impossible for a single galaxy 
this was unlikely for an ensemble of lenses. In a subsequent analysis, Dobke 
and King [55] did this the other way round; they took the value of Hq = 72 ± 
8 km s _1 Mpc -1 in [45] and deduced that the overall mass slope index in time- 
delay lenses had to be 0.2-0.3 steeper than isothermal. If true, this is worrying 
because investigation of a large sample of SLACS lenses with well-determined 
mass slopes [52] reveals an average slope which is nearly isothermal. 

More recent analyses, including information available since then, may be 
indicating that the lens galaxies in some earlier analyses may indeed, by bad 

bananas. 
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luck, be unusually steep. For example, Treu and Koopmans [159j find that 
PG1115+080 has an unusually steep index (0.35 steeper than isothermal) yield- 
ing a > 20% underestimate of Hq. The exact value of Hq from all eighteen 
lenses together is a rather subjective undertaking as it depends on one's view 
of the systematics in each lens and the degree to which they have been re- 
duced to random errors. My estimate on the most optimistic assumptions is 
66 ±3 km s _1 Mpc -1 , although you really don't want to know how I got thij^l. 

A more sophisticated meta-analysis has recently been done in [102] using a 
Monte Carlo method to account for quantities such as the presence of clusters 
around the main lens galaxy and the variation in profile slopes. He obtains 
(68 ± 6 ± 8) km s _1 Mpc -1 . It is, however, still an uncomfortable fact that 
the individual Hq determinations have a greater spread than would be expected 
on the basis of the formal errors. Meanwhile, an arguably more realistic ap- 
proach [127j is to simultaneously model ten of the eighteen time-delay lenses 
using fully non-parametric models. This should account more or less automat- 
ically for many of the systematics associated with the possible galaxy mass 
models, although it does not help us with (or excuse us from further work to 
determine) the presence of mass sheets and their associated degeneracies. The 
result obtained is 72+* 2 km s -1 Mpc -1 . These ten lenses give generally higher 
H values from parametric models than the ensemble of 18 known lenses with 
time delays; the analysis of these ten according to the parametric prescriptions 
in Appendix lAl gives Hq = 68.5 rather than 66.2 km s _1 Mpc -1 . 

To conclude; after a slow start, lensing is beginning to make a useful contri- 
bution to determination of Hq, although the believable error bars are probably 
similar to those of local or CMB methods about eight to ten years ago. The 
results may just be beginning to create a tension with other methods, in the 
sense that Hq values in the mid-sixties are preferred if lens galaxies are more or 
less isothermal (see |76j for discussion of this point). Further work is urgently 
needed to turn systematic errors into random ones by investigating stellar dy- 
namics and the neighbourhoods of galaxies in lens systems, and to reduce the 
random errors by increasing the sample of time delay lenses. It is likely, at 
the current rate of progress, that < 5% determinations will be achieved within 
the next few years, both from work on existing lenses and from new measure- 
ments of time delays. It is also worth pointing out that lensing time delays 
give a secure upper limit on Hq, because most of the systematic effects associ- 
ated with neglect of undetected masses cause overestimates of Hq; from existing 
studies Hq > 80 km s _1 Mpc -1 is pretty much ruled out. This systematic of 
course makes any overall estimates of Hq in the mid-sixties from lensing very 
interesting. 

One potentially very clean way to break mass model degeneracies is to dis- 
cover a lensed type la supernova [10311104] . The reason is that, as we have seen, 
the intrinsic brightness of SNe la can be determined from their lightcurve, and it 
can be shown that the resulting absolute magnification of the images can then 
be used to bypass the effective degeneracy between the Hubble constant and 

17 Full details are in Appendix El 
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the radial mass slope. Oguri et al. [104) and also Bolton and Buries [13] discuss 
prospects for finding such objects; future surveys with the Large Synoptic Sur- 
vey Telescope (LSST) and the SNAP supernova satellite are likely to uncover 
significant numbers of such events. With modest investment in investigation 
of the fields of these objects, a < 5% determination of Hq should be possible 
relatively quickly. 

4.2 The Sunyaev— Zel'dovich effect 

The basic principle of the Sunyaev-Zel'dovich (S-Z) method |147j . including its 
use to determine the Hubble constant [141] , is reviewed in detail in [T21 12U] . It 
is based on the physics of hot (10 8 K) gas in clusters, which emits X-rays by 
bremsstrahlung emission with a surface brightness given by the equation 

^MiW/^ a (14) 

(see e.g. [H]), where n e is the electron density and A e the spectral emissivity, 
which depends on the electron temperature. 

At the same time, the electrons of the hot gas in the cluster Compton up- 
scatter photons from the CMB radiation. At radio frequencies below the peak 
of the Planck distribution, this causes a "hole" in radio emission as photons 
are removed from this spectral region and turned into higher-frequency photons 
(see Figure [5]). The decrement is given by an optical-depth equation, 

=I J n e a T y(v,T e )dl, (15) 

involving many of the same parameters and a function which depends on 
frequency and electron temperature. It follows that, if both &x and AI(x) can 
be measured, we have two equations for the variables n e and the integrated 
length in through the cluster and can calculate both quantities. Finally, if we 
assume that the projected size l±_ of the cluster on the sky is equal to In, we 
can then derive an angular diameter distance if we know the angular size of the 
cluster. The Hubble constant is then easy to calculate, given the redshift of the 
cluster. 

Although in principle a clean, single-step method, in practice there are a 
number of possible difficulties. Firstly, the method involves two measurements, 
each with a list of possible errors. The X-ray determination carries a cali- 
bration uncertainty and an uncertainty due to absorption by neutral hydrogen 
along the line of sight. The radio observation, as well as the calibration, is 
subject to possible errors due to subtraction of radio sources within the cluster 
which are unrelated to the S-Z effect. Next, and probably most importantly, 
are the errors associated with the cluster modelling. In order to extract pa- 
rameters such as electron temperature, we need to model the physics of the 
X-ray cluster. This is not as difficult as it sounds, because X-ray spectral in- 
formation is usually available, and line ratio measurements give diagnostics of 
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Figure 8: S-Z decrement observation of Abell 697 with the Ryle telescope in 
contours superimposed on the ROSAT grey-scale image. Reproduced from \69}/ . 
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Table 2: Some recent measurements of Hq using the S-Z effect. Model types are 
f3 for the assumption of a fi-model and H for a hydrostatic eguilibrium model. 
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physical parameters. For this modelling the cluster is usually assumed to be in 
hydrostatic equilibrium, or a "beta-model" (a dependence of electron density 
with radius of the form n(r) = no(l + r 2 /r 2 )~ 3 ^/ 2 ) is assumed. Several recent 
works |137i 15J relax this assumption, instead constraining the profile of the 
cluster with available X-ray information, and the dependence of Hq on these 
details is often reassuringly small (< 10%). Finally, the cluster selection can be 
done carefully to avoid looking at cigar-shaped clusters along the long axis (for 
which l± 7^ iii) and therefore seeing more X-rays than one would predict. This 
can be done by avoiding clusters close to the flux limit of X-ray flux-limited 
samples, Reese et al. |121] estimate an overall random error budget of 20-30% 
for individual clusters. As in the case of gravitational lenses, the problem then 
becomes the relatively trivial one of making more measurements, provided there 
are no unforeseen systematics. 

The cluster samples of the most recent S-Z determinations (see Table [2]) 
are not independent in that different authors often observe the same clusters. 
The most recent work, that in [T5j is larger than the others and gives a higher 
Hq. It is worth noting, however, that if we draw subsamples from this work and 
compare the results with the other S-Z work, the H values from the subsamples 
are consistent. For example, the H derived from the data in [15] and modelling 
of the five clusters also considered in is actually lower than the value of 
66 km s _1 Mpc -1 in [55] , 

It therefore seems as though S-Z determinations of the Hubble constant are 
beginning to converge to a value of around 70 km s _1 Mpc -1 , although the 
errors are still large and values in the low to mid-sixties are still consistent with 
the data. Even more than in the case of gravitational lenses, measurements of 
Ho from individual clusters are occasionally discrepant by factors of nearly two 
in either direction, and it would probably teach us interesting astrophysics to 
investigate these cases further. 

4.3 Gravitational waves 

One topic that may merit more than a brief paragraph in a few years' time is 
the study of cosmology using gravitational waves. In particular, a coalescing 
binary system consisting of two neutron stars produces gravitational waves, and 
under those circumstances the measurement of the amplitude and frequency of 
the waves determines the distance to the object independently of the stellar 
masses |138j . This was studied in more detail in [22] and extended to more 
massive black hole systems 58, 29] . More massive coalescing signals produce 
lower-frequency gravitational wave signals which can be detected with the pro- 
posed LISA space-based interferometer. The major difficulty is obtaining the 
rcdshift measurement to go with the distance estimate, since the galaxy in which 
the coalescence event has taken place must be identified. Given this, however, 
the precision of the Hq measurement is limited only by weak gravitational lcns- 
ing along the line of sight, and even this is reducible by observations of multiple 
systems or detailed investigations of matter along the line of sight. Ho deter- 
minations to ~ 2% should be possible in this way after the launch of the LISA 
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satellite in 2015. 



5 Conclusion 



The personal preference of the author, if forced to calculate a luminosity from a 
flux or a distance from an angular size, is to use the value of Hq ~ 70 - 73 km s _1 Mpc 
which is fast becoming standard. No one determination is individually com- 
pelling, but the coincidence of the most plausible cosmological determinations 
and some local determinations, together with the slow convergence of the one- 
step methods of gravitational lensing and the Sunyaev-Zel'dovich effect, is be- 
ginning to constitute a consensus that it is hard to visualise being overturned. 
Most astronomers, if forced to buy shares in a particular value, would be likely 
to put money into a similar range. It is, however, possibly unwise to sell the 
hedge fund on a 10% lower value just yet, and it is important to continue the 
process of nailing the value down - if only for the added value that this gives in 
constraints that can then be applied to other interesting cosmological parame- 
ters. Barring major surprises, the situation should become very much clearer in 
the next five to ten years. 
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A A Highly Subjective Summary of the Grav- 
itational Lens Time Delays to Date and the 
Derived Values of the Hubble Constant 

All discussion is based on parametric models (though see the text for non- 
parametric models). For B0218+357 the mass-model is reasonably well known, 
thanks to an Einstein ring and VLBI structure |176| but the galaxy centre as- 
trometry is not. York et al. [178] find 70±5 (2a) based on an ACS determination 
of the galaxy centre, and Wucknitz et al. [176j find 78 ± 6 (2a) using the ring 
and VLBI structure. I double the errors and take 74 ± 5. For HE0435-1223 the 
image positions of [174] with the lens redshift of [HH] and the longest time de- 
lay of [79] , together with an isothermal sphere model with external shear, gives 
Hq = 90, similar to the value inferred by Kochanek et al. A 20% error for the 
lack of knowledge of the mass profile and another 20% for the presence of nearby 
galaxies are added in quadrature to the time delay error to give Hq = 90 ± 26. 
For SBS0909+532 the system is very underconstrained, especially as the flux 
ratios are not well known, possibly due to extinction |166| . Lehar et al. [55J 
find a predicted time delay of 49 ± 38 days for Hq — 70, by considering the 
allowed space where x 2 < 1, the time delay of 45 days |166j translates into 
Hq = 76 ± 59; one might argue that the errors should be slightly smaller but 
this lens is not in any case going to contribute significantly to the overall to- 
tal. In RXJ0911+0551 the situation is complicated by the presence of a nearby 
cluster. Hjorth et al. [SHj obtain 71 ± 9 (2a) based on a velocity dispersion for 
the cluster from [74] (but see also discussion in [77]). As with B0218+357 I 
double the errors and take 71 ±9. In FBQ0951+2635 an isothermal model gives 
60^^2 EE] ■ The same source claims that in this lens the variation of Hq with 
mass slope is relatively weak, but I nevertheless add an extra 20% to the error 
and take 60 ± 15. For Q0957+561, a system with a huge cluster contribution, 
probably the most compelling analysis (of the vast literature on this object) 
is that in 70 where the lensed host galaxy is used to disentangle the cluster 
shear from the internal shear produced by the lensed galaxy shape. Assuming a 
cluster convergence of 0.2, using the range of Keeton et al.'s sample models and 
increasing the errors by 20% gives 90 ± 25. For SDSS J1004+4112 the cluster 
contribution makes modelling so problematic that it is difficult to sort out sensi- 
ble Hq estimates from simple models. However, once the cluster mass profile is 
thoroughly understood using the large number of constraints and the remaining 
time delays are known, this may give the best Hq estimate of all [12] ■ HE 11 04- 
185 is also a bizarre system in which the closer image is brighter; Lehar et al. [55] 
predict much larger time delays than are actually found, and typical isothermal 
models plus the time delay in |101j give Hq ^95. For consistency I adopt this 
value, together with a 20% error added in quadrature to the time delay error to 
reflect the fact that we do not know the galaxy mass profile, giving 95 ± 25. In 
PG1 115+080 traditional fits of isothermal models give Hq in the forties |135j . 
but a more recent stellar dynamical measurement .159] suggests a steeper-than- 
isothermal profile and with the time delay in [135j gives Hq = Wt) 2 ^, so I 
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adopt 59 ± 12. In RXJ1131-1231 almost any smooth isothermal model fails 
(and H < 10), so Morgan et al. [98; suggest that a large piece of substructure 
or satellite galaxy falls close to one of the images. While plausible, this means 
that the system is almost impossible to model for Ho- There is a measured time 
delay for JVAS B1422+231 [110) but in my view it needs confirmation before 
being used for H determination. SBS 1520+530 is modelled in [IS], who also 
provide the time delay, by an isothermal model together with consideration of a 
nearby cluster. I add 20% to their error, due to the isothermal assumption, and 
get 51 ± 13. In B1600+434 an isothermal model gives 60±^ (random errors ol , 
but adjusted for a (Q m ,flA — (0.3,0.7) model). I add 20% for systematics in- 
volved in lack of knowledge of the mass profile and get 60 ± 20. In B1608+656 
Koopmans et al. [83] assemble an impressive array of data, including three time 
delays, stellar velocity measurements and Einstein ring fitting to get 75lg, to 
which I add another 10% due to a feeling of unease about the lens being two 
interacting galaxies and consequent effects on the mass profile, and take 75 ±10. 
For SDSS J1650+4251 I adopt the value of H = 52 from the time delay and 
isothermal models of |167j . plus an extra 20% error for the mass slope, and use 
H = 52 ± 11. In PKS1830-211 models in [173] using their improved galaxy 
centre position give 44 ±9, to which I add the now traditional 20% for ignorance 
of the mass slope and 10% (cf. B1608+656) for the possibility of associated sec- 
ondaries [25], to give 44 ± 13. Finally, for HE2149-2745, a simple SIE model 
gives H — 49, in agreement with the average of 48^ quoted in [73] on anal- 
ysis of this plus four other lenses; Burud et al. [17] derive a higher value due 
to restrictions on assumed ellipticity of the lens galaxy. Adding the usual 20% 
for ignorance of the mass slope to the time delay error we obtain 49 ± 11. At 
last, the overall average is 66 ± 3 (or 61 ± 4 without B0218+357). We should 
also consider the effect of mass along the line of sight, which is likely to add a 
further 5-10% to the error budget of each lens system and which will probably 
not reduce as ^/n; much of it is systematic if the error is normally in the sense 
of ignoring nearby groups of matter which add to the local convergence. 

It should be emphasised that this is an attempt to work out one observer's 
view of the situation based on parametric models, and that the reader will have 
noticed a lot of subjective judgements and arbitrary manipulation of errors going 
on, but you did ask. 
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